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Abstract 

Space-filling assemblies of athermal hydrophobic particles floating at an air- water interface, called 
particle rafts, are shown to undergo an unusual phase transition between two i.e., a low density 
'less-rigid and a high density 'more-rigid' amorphous states as a function of particulate number 
density ($>). The former is shown to be a capillary-bridged solid and the later a frictionally-coupled 
one. Simultaneous studies involving direct imaging as well as measuring its mechanical response to 
longitudinal and shear stresses show that the transition is marked by a subtle structural anomaly 
and a weakening of the shear response. The structural anomaly is identified from the variation 
of the mean coordination number, mean area of the Voronoi cells and the spatial profile of the 
displacement field with <£. The weakened shear response is related to local plastic instabilities 
caused by the depinning of the contact-line of the underlying fluid on the rough surfaces of the 
particles. 

PACS numbers: 46.55.+d 



I. INTRODUCTION 



In crystalline solids in thermodynamic equilibrium, the onset of rigidity is a direct con- 
sequence of the appearance of long range positional order and the broken continuous trans- 
lational symmetry. No such overarching principles are known to govern amorphous solids 
which are frequently out of equilibrium. For example, the onset of rigidity in athermal gran- 
ular systems has been a topic of great recent interest pQ . Such studies explore the transition 
of a system from a state of zero rigidity (unjammed state) to that of a finite one (jammed 
state) [2H3]. But only a few have explored further phase transitions that may exist within 
the jammed state of a system [B]. In this paper we study the mechanical response of a rigid 
but amorphous particle raft to compressive (longitudinal) and shear (transverse) stresses for 
varying particulate number density (<&). These rafts are space-filling assemblies of athermal 
hydrophobic particles floating at an air-water interface and have properties common to elas- 
tic [7] and granular [8j solids. Its mechanical response and structural reorganization reveal 
anomalies that are suggestive of a phase transition between two amorphous states, i.e., a 
low-density 'less-rigid' state and a high-density 'more-rigid' state. 

II. EXPERIMENTAL DETAILS 

A. Particle raft: A rigid and athermal model system 

Surface tension assists hydrophobic particles which are denser than water to float on 
it. They deform the otherwise flat liquid surface [9], [10] under gravity, thereby generating 
long-range inter-particle attraction and form particulate-clusters (TTJ [T2]. The individual 
clusters show solid like properties, i.e., they retain their shape and have a finite rigid- 
ity. However, for small areal coverage, these clusters are sparsely distributed. Hence, at 
length scales comparable to the system size the collective mechanical property of the floating 
clusters is governed by the intervening liquid. Upon compression, i.e., increasing the areal 
coverage, the individual clusters fuse to form a contiguous system spanning structure. The 
short-range inter-particle interactions depend on the surface roughness and wettability of 
the particles and are either attractive or repulsive [TSJ E33]- The poly-dispersity and the 
athermal nature of the particles make the raft amorphous. 
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B. Experimental protocol to prepare the particle raft and its structural charac- 
terisation 

The following experimental protocol is used to prepare the space filling structure of the 
raft. 

(i) The hydrophobic (coated with FluoroPel PFC M1104V-FS from Cytonix LLC) silica par- 
ticles of average diameter (2a) are initially sprinkled on the air-water interface in a Langmuir 
trough. These particles form disjoint particulate clusters (see FigjlTa)). This state of the 
system is defined as a ' patchy' state. For data presented in this paper, unless otherwise 
mentioned, 2a =0.5mm (polydispersity is 15% and density of silica is 2500 Kg/m 3 ). 

(ii) These clusters are then brought within the interaction distance, i.e., capillary length 
(L c ~ 2mm), [HI Q3] by moving the motorized Teflon barriers of the trough inwards. This 
constitutes the first compression cycle (cl). As a result of this compression the disjoint clus- 
ters coalesce to form a system-spanning quasi-two dimensional raft [7J El H3] • The inward 
motion of the Teflon barriers is stopped just before the out-of-plane deformations (wrinkling) 
[7] of the raft appear. The resulting homogeneous 'compressed' state of the particle raft is 
shown in Fig. [I](b). 

(iii) The barriers are then moved outward in the first expansion cycle (el) until they detach 
from the raft (see movie 'MSl.avi' in Supplementary information). Figure [l|c) shows the 
micrograph of the resulting relaxed, yet rigid, state of the raft. Further cycles of com- 
pression and expansion transform the system between the 'compressed' and the 'relaxed' 
states. 

The number density of the particles (<£>) between the barriers defined as, $ = 
Nira 2 /(L x L y ), where L x and L y (= 140mm) are the length and width of the raft, re- 
spectively is chosen to be the relevant control parameter. This is guided by the literature 
on jamming transition [H [2]. The reported measurements are made in a range below the 
density where folds, i.e., out of plane distortions [7] occur (see movie 'MS2.avi' in Suppl. 
Information). The lack of distinct features beyond the third peak in the radial density pair- 
correlation function, g(r), in Fig. []Jd), illustrates that the raft remains amorphous in 
both the 'compressed' (see Fig. [l](b)) and the 'relaxed' (see Fig. [l](c)) states. 
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C. Experimental setup 



The schematic of the experimental set-up is shown in Fig. [2] The system is 'compressed' 
(or expanded) by moving two motorized Teflon barriers along x in steps of 0.05mm. The 
system is sheared sinusoidally in the x-direction with a hydrophobized microscope cover- 
slip attached to a piezo-stage (PI-517.3CL) with a amplitude (w°)=0.05mm and a frequency 
(f)=20Hz. A stainless-steel cantilever and a parallel plate capacitor is used to measure shear 
(u xy ) and longitudinal (a xx ) stresses, respectively. Details of the experimental methods and 
the parameters used for calculating stresses are described in Appendix A. Additionally the 
system is imaged under no shear for each position of the barriers. These images are then 
analyzed to obtain the mean coordination number (Z) and to generate the Voronoi diagram 
from which mean cell area (A) is calculated. 

D. Measurement of mechanical response 

The mechanical response of this system to external stresses (longitudinal and shear) is 
described in this paper in terms of a spatially averaged effective longitudinal and shear 
moduli (Ka an d G, respectively) that are defined as: Ka = 8<r xx /5u xx and G = <J xy /u xy , 
where u xy = u° x /D is the shear strain, D = 60mm is the distance of the cantilever from 
the shear-launching microscope cover-slip and Su xx = —5L X /L X = 8<&/<&o is the incremental 
compressive strain, where L x is the length of the raft in the 'relaxed' ('compressed') state 
for an compression (expansion) run and $o is the density corresponding to L x . Both a xx 
and a xy are measured with an instrumental resolution of 1 mPa. The noise observed in the 
data is intrinsic to the system and is a signature of finite size effect of the system. Thus the 
longitudinal modulus (Ka) is calculated by numerically differentiating smoothed (using a 
polynomial fit) a xx with respect to $. We note that the elastic moduli are used here as an 
intuitive measure of the stress-transmission defined above [15] . 

Figure [3] (a) and (b) show the variation of a xx and a xy with $ for cl (open symbols) 
and el (filled symbols). The first compression cycle (cl) starts from the patchy state (see 
Fig. [T] (a)) where both a xx and a xy are zero, while the expansion cycle (el) starts from 
the 'compressed' state (see Fig. [T| (c)). During the first compression, system-size spanning 
stress-bearing networks form around $ = 0.74, marked by rapidly growing longitudinal and 
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shear stresses. However, during the subsequent expansion cycle towards the 'relaxed' state 
(<3> ~ 0.69, as in Fig. 1(c)), both a xx and a xy remain finite. 

For the subsequent cycles (the first expansion onwards) the observed variation of the 
stresses separates into two regions: Region I for $ < 0.74, where a xy decreases and region 
II for $ > 0.74, where a xy increases rapidly with $, whereas a xx increases monotonically and 
non-linearly with <3> in both regions. The variations in the stress transmission characteristics 
are illustrated through the computed Ka and G, shown in Fig. [I] (a) and (b), respectively. 
The ratio, Ka/G, shows a pronounced but inhomogeneously jagged cusp around $ = 0.74 
(see Fig. [4j^c)) [16J. The transition region is broad and shown as a shaded region in the 
figures. 

E. Structural rearrangement 

The displacement field associated with the particles in response to the barrier movement is 
calculated by digitally cross-correlating an image corresponding to a specific $ to a later one 
obtained after displacing the barriers by 0.25mm. The details of computing the displacement 
field is given in |T7|. The right panel of Fig. [5] shows the typical displacement field during 
the first expansion cycle. The field is inhomogeneous and exhibits a large number of complex 
cellular features for $ < 0.74. This is reflected in the broad distribution of the probability 
distribution function, P{0), of the argument (8) associated with the displacement vectors 
(s), i.e., 6 = tan~ 1 (s y /s x ) [18] shown in the left panel of Fig. [5] However, for $ > 0.74 
the field becomes homogeneous and P(0) narrows considerably. The spatial resolution of 
the above imaging method is 0.05 mm and it limits us from calculating the displacement 
field associated with the shear deformation. 

The inhomogeneity of the displacement field at low densities ($ < 0.74) arise from 
mechanical instabilities which can be interpreted as plastic events at the scale of the particle 
size (see movie 'MS3.avi' and Fig. [6]). In order to study the plasticity associated with 
the motion of the particles, each of them are individually tracked. Representative particle 
displacements (Aj) measured from the 'relaxed' state is plotted as a function of a xx in Fig. [6] 
(a) for the second compression (c2) cycle. The jaggedness of the graph, i.e., sudden jumps in 
Aj's, is associated with plasticity. A measure of the total number of plastic events associated 
with the stress cycling is obtained from the hysteresis of the cumulative strain parameter 
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(A),= J^Ai/L°, where the summation is carried out over all the particles. Figure [6^b) 
shows the variation of A with a xx for successive compression (open symbols) and expansion 
(filled symbols) cycles. 

The images of the raft at various $'s are used to further explore possible structural 
signatures of the transition. Figure [7] (a) shows an increase in the mean coordination number 
Z (= N J^ +e 2irlg(l)dl), where e is chosen to match the first peak in g{r) as a function 
of $. The corresponding variation of the Voronoi cell area A normalized with the particle 
area (vra 2 ), an alternative measure of the local structure and compactivity [TH], is plotted in 
Fig. [7] (b). Both Z and < A > /rca 2 varies linearly with $ in region I, but however start 
to deviate from the initial slope (shown by the extrapolated dotted line) in the transition 
region. It suggests that the measured weakening of the shear modulus is accompanied by 
structural reorganizations at the length scale of nearest neighbors throughout the system. 



III. DISCUSSION OF RESULTS 



The choice of an order parameter in a disordered system is usually not unique and is often 
subjective [20]. In our experiments we find the mean coordination number Z is a sensitive 
parameter to describe the transition. The dependence of Ka and G on Z is shown in Fig. 
[8] (a) and (b) respectively for el and c2. The system exhibits two terminal states, (i) 'less- 
rigid' state (Ka ~ 20Pa and G ~ 10 Pa) and (ii) a 'more-rigid' state (Ka ~ 2000Pa and 
G ~ 200Pa). The transition from one state to the other is marked by decrease in G around 
Z ~ 3.6, although, Ka changes monotonically with increasing Z. Hence, the anomaly 
observed is in the variation of the shear compliance, i.e., 1/G, reminiscent of the behavior 
of the magnetic susceptibility at a spin glass transition in disordered magnets |21j . 

The cusp in Ka/G (see Fig. [4] (c)) is a robust feature of this transition. This is a 
signature of the difference in the density dependence of Ka and G, reminiscent of, but 
distinct from, the power-law divergence observed in simulations of the jamming transition 
in frictionless granular materials [1] . However, one must emphasize that the phenomenology 
observed here is far from the original definition of athermal jamming for which only hard 
core interaction is considered and where the transition is between rigid and non-rigid states. 
In the conventional sense the 'relaxed' state of the system is already in a jammed state. The 
observations should be viewed as a transition between two jammed states whose differences 
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are discussed below. 

A. Microscopic mechanism associated with the transition 

In order to understand the microscopic mechanism associated with the transition, we 
investigated the dependence of G on the particle size (see inset of Fig. [8] (b)). For the 
'less-rigid' state (filled triangles), G oc a -1 . This is consistent with a capillary bridging 
mechanism (G ~ 7/a, where 7 is the surface tension of the liquid) through the pinned 
contact lines of the liquid on the particles [13]. In the 'more-rigid' state the contact 
friction dominates the inter-particle interaction. The measured shear modulus (G) scales 
with the number density of contacts, which is proportional to the number density of par- 
ticles, i.e., G oc N/(L x L y ) ~ $/a 2 . It assumes that the roughness scale is independent of 
the particle size. Since this happens over a narrow range of $, one obtains G oc a~ 2 as in 
conventional elasticity in a quasi-2D system [22]. We hence conclude that regions I and II 
represent a capillary-bridged solid and a frictional solid, respectively. 

B. Time dependent effects and stress cycling in the particle raft 

The structural relaxation time of the particle raft, r = 12irr]a 3 /mv 2 ~ 1000s, where 
m = 0.2fig and v = 50/xm/s are the mass and velocity of a single particle respectively and 
77 = ImPas is the viscosity of the liquid. In calculating the relaxation time, we assume 
that the effective temperature of the system (T e ff) is related to the kinetic energy of the 
particle, i.e., T e ff = mv 2 /2Kg, where Kb is the Boltzmann's constant. The system thus 
takes a long time to reach its equilibrium state. This would mean that the system would 
typically find itself in a kinetically arrested state and would show a dependence on the 
history of the paths via which a given state of the system is reached. Indeed, measurements 
above 1 mHz probes a frequency-dependent rigidity of the material as shown in Fig. [9]^a) for 
$ = 0.77. Moreover, creep effects (time-dependence) in a xy are more pronounced in region 
II than in region I, as shown in Fig. [9] (b). Interestingly, the stress allows the 'more-rigid' 
frictional solid to explore various metastable states and hence the system shows a marked 
'creep', i.e. a temporal variation, in <j xy . This is absent for the 'less-rigid' capillary solid 
suggesting that it is closer to a deeper metastable minimum. Strong history dependence is 
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also seen in the variation of a xx and a xy measured for subsequent compression and expansion 
cycles which broadly follow the trend (see fig. |3]) but with greatly reduced hysteresis (see 



Fig. 10, the data shown in the figure is for particles whose diameter is 1mm), analogous to 
residual densification observed in amorphous materials |23j . 



IV. CONCLUSION 



We have shown that a rigid particle raft can undergo a phase transition from a 'less-rigid' 
low density state to a 'more-rigid' high density state as a function of particulate number 
density. The transition is marked by a weakening of the shear modulus. The measured 
shear modulus which distinguishes the two states of the system, i.e., 'less-rigid' and 'more- 
rigid', arises from the restoring force of pinned contact lines in the former case and from 
particle contacts through a frictional coupling in the latter. The weakening of the shear 
modulus observed in the crossover region is thus attributable to a reduction of restoring 
force arising from plastic events caused by the mechanical instabilities and the associated 
depinning of the contact lines [2H [25] . 

Although the results presented in the present paper are specific to particle rafts, the weak- 
ening of the shear response and local reorganization seen in the present experiment, have 
also been observed in network glasses [251 EZ]- We hope that the present experiment will 
provide some insight into the complex phenomenon of pressure induced phase transitions in 
amorphous solids (23J, 1281432] . 

The authors thank M. Bandi, M. Cates, N. Menon, L. Mahadevan, S. R. Nagel and A. 
Yethiraj for helpful discussions. 



Appendix A: Experimental set-up 

The schematic of the experimental set-up is shown in Fig. [2] Hydrophobic silica particles 
are sprinkled on the air-water interface in a Langmuir trough (300mm x 145mm x 25mm) 
made of Teflon. These particles coalesce and form disjoint (patchy) clusters. The two 
motorized Teflon barriers of the trough are used to fuse these clusters into a homogeneous 
'compressed' or 'relaxed' state by moving the barriers inward or outward in x-direction 
respectively. The system is illuminated uniformly from the bottom using an LCD monitor 
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through a glass window sealed in the trough and a camera above records the image in the (x— 
y) plane at each barrier step (= 0.05mm). The system is sheared in the x-direction by moving 
a hydrophobic cover-slip (60mm x 24mm x 0.15mm) sinusoidally u x = u x cos(2irvt), where 
u° x = 0.05mm and frequency v = 20Hz, connected to a three-axis piezo-stage (PI-517.3CL). 
A stainless steel cantilever is kept at D — 60mm in the y-direction from the shear launching 
cover-slip as marked by (i) and shown at the left-bottom of the Fig. [2] with dimensions 
mentioned. The hydrophobic cover-slips attached to piezo stage and the cantilever makes 
contact with the particles at the air-water interface. The shear stress {(T xy ) is measured from 
the lateral displacement of the cantilever. The displacement is sensed by a fiber optic-based 
displacement sensor (MTI-2000 FOTONIC) from the interference of light that is reflected 
from the mirror glued on the cantilever at a distance of 25mm from the fixed support. The 
output of the displacement sensor is fed to a lock-in amplifier (SR830). A parallel plate 
capacitor, used as a force-sensor, is marked by (ii) and also shown with dimensions in 
the right-bottom of Fig. [2] It is made of 0.2 mm thick flexible polymeric sheet and kept 
at a distance of 30mm from the cantilever in the ^-direction. The metal electrodes are 
glued at the bottom of the two plates of the force-sensor with a separation of 3mm and 
are immersed in water. A capacitance-bridge (16 15- A General Radio) is used to detect the 
change in the plate separation during compression (or expansion) from the lateral deflection 
of the individual capacitor plates. The output of the capacitance-bridge is fed to another 
lock-in amplifier (SR830) which drives the bridge at lOOKHz. The stresses, shear (<J xy ) 
and longitudinal (a xx ), are calculated using "bending-of-a-beam" formula using lateral 
displacement of the cantilever in the former case and the plate deflection in the latter. The 
formula used in general is given by: the stress a = QREI /xq(3/o — xo)a |33j where R is the 
lateral displacement, E is Young's modulus, I = bb? j\2 is the second moment of area (b 
and h are width and thickness of the detecting object, respectively), l is the length of the 
object, Xq is point of detection from the fixed end and a is an effective area of contact with 
the particles. All these parameters used in the calculation of stress' are tabulated in table [j] 
for the cantilever and the parallel plate force sensor. 
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length (In) 
(mm) 


width (b) 
(mm) 


thickness (h) 
(mm) 


Young's 
modulus (E) 
(GPa) 


point of de- 
tection (xo) 
(mm) 


effective area 
(a) 

(mm 2 ) 


Cantilever 


90 


13 


0.13 


200 


25 


30 


Parallel plate 
force-sensor 


60 


25.4 


0.2 


2 


60 


12.7 



TABLE I: Parameters used in the calculation of a xx and a xy for the cantilever and the parallel 
plate force sensor 
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(a) Patchy State 
■ 




(b) Compressed State 




FIG. 1: Micrographs of hydrophobized silica particles at the air- water interface in the form of (a) 
disjoint clusters (patchy state) and (b) a homogeneous 'compressed' state of the raft in presence and 
(c) absence ('relaxed' state) of longitudinal stress. The scale bars, 10mm. (d) The left (right) panel 
shows the plot of the radial density pair correlation function, g(r), of the 'relaxed' ('compressed') 
state of the system. 
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FIG. 2: Experimental setup: hydrophobized microscope cover-slip attached to a piezo-stage oscil- 
lates along x producing a sinusoidal shear deformation. The cantilever with an attached microscope 
cover-slip, placed along y at a distance (Z))=60mm away from the piezo-stage measures the shear 
stress (cr xy ). A parallel-plate capacitor, marked as the force-sensor, is used to measure the 
longitudinal stress (a xx ). 
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FIG. 3: Variation of longitudinal stress (a xx ) and shear stress (cr xy ) are plotted as a function of 
<3? in (a) and (b) respectively for cl (open symbols) and el (filled symbols) cycles. The arrows 
in (a) shows the direction in which <I> changes. The region over which u xy shows a softening 
is marked as region I. In region II, a xy increases rapidly. The different states of the system 
i.e., patchy, 'compressed' and 'relaxed' are marked in the figure. The transition region separating 
I and II is shaded. 
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FIG. 4: The variation of (a) longitudinal modulus (Ka), (b) shear modulus (G), (c) and the ratio 
Ka/G as a function of for the first expansion cycle. The longitudinal modulus (Ka) is calculated 
by numerically differentiating smoothed (using a polynomial fit) a xx with respect to 
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FIG. 5: Displacement field of the particles and the probability distribution of angles, P{&) as- 
sociated with the displacement vector (s) measured during the first expansion cycle is shown in 
right and left panel respectively. Here 9 = ian _1 (^-) For values of <3? < 0.74 the displacement 
field is inhomogeneous and exhibits a large number of complex cellular features (right panel) . This 
is reflected in the broad distribution of P{9) (left panel). However, for values of > 0.74, the 
displacement field becomes spatially correlated over large distances and P{0) narrows considerably. 
The arrow indicates the progressive sequence of $'s, i.e, the system moves from a state of high 
density, i.e., 3> = 0.77 to a state of low density, i.e., = 0.69. 
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FIG. 6: (a)Typical displacement curves of a few particles (Aj), measured from its position 
in the 'relaxed' state as a function of a xx . For visual clarity the curves are shifted vertically. 
The corresponding scale bar of the displacement is shown in the figure. The data is plotted for a 
compression cycle, (b) The variation of the cumulative strain ( A = ^ Ai/L x , where the summation 
is over all the particles) with a xx for successive compression (open symbols) and expansion (filled 
symbols) cycles. The uncertainty in the measurement of Aj and A is shown by the accompanying 
band. 
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0.68 0.70 0.72 0.74 0.76 0.78 

FIG. 7: The variation of (a) the mean coordination number Z and (b) the mean Voronoi cell area 
A normalized with the particle's area (vra 2 ) as a function of The solid lines show their linear 
variation in region I. The dotted lines which are a linear extrapolation of initial variation highlights 
the later deviation of the graphs. The data presented here is for el. 
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FIG. 8: Variation of (a) Ka and (b) G as a function of Z. The filled symbols correspond to el 

while the open ones correspond to c2. The shaded regions indicate the spread in the data. The 

'less-rigid' (schematically shown as a capillary bridged interaction in (a)) and 'more-rigid' states 

(schematically shown as a friction dominated interaction in (a)) are marked in the figure. The 

inverted cusp in G versus Z corresponds to the slipping of the contact line (schematically shown in 

(b) with a greater width of a sliding contact line). The inset of (b) shows the variation of G of the 

'less-rigid' (triangles) and 'more-rigid' (square) states of the raft as a function of particle diameter 

(2a). The solid lines passing through the data show a -1 and a -2 variation of G for the 'less-rigid' 
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and 'more-rigid' states, respectively. 
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FIG. 9: (a) Frequency dependence of G measured at <!> = 0.77. (b) The time dependence (creep 
effect)of a xy is shown for different <5's. 
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FIG. 10: The figure shows the variation of a xx axi&<j xy for various compression (open symbols) 
and expansion (filled symbols) cycles as a function of for particles of diameter 1 mm. The 
sequence of expansion and compression are as follows, first compression (data not shown for 
reasons discussed in the main text) : first expansion (el: filled triangles), second compression 
(c2: open triangles), second expansion (e2: filled circles), third compression (c3: open circles), 
third expansion (e3: filled stars), fourth compression (c4: open stars). The repeated compression- 
expansion cycles result in annealing of the raft and thus later cycles show a reduced but finite 
hysteresis. 
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